Instant preheating mechanism and UHECR 
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Top-down models assume that the still unexplained Ultra High Energy Cosmic Rays (UHECR's) 
are the decay products of superheavy particles. Such particles may have been produced by one of the 
post-inflationary reheating mechanisms and may account for a fraction of the cold dark matter. In 
this paper, we assess the phenomenological applicability of the simplest instant preheating framework 
not to describe a reheating process, but as a mechanism to generate relic supermassive particles as 
possible sources of UHECR's. We use cosmic ray flux and cold dark matter observational data to 
constrain the parameters of the model. 
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I. INTRODUCTION 

The possible observation of Ultra High Energy Cos- 
mic Ray (UHECR's) events with primary energies above 
10 20 eV [1] constitute one of the most intriguing puzzles 
in astroparticle physics (see, for example, [2]), although 
their origin and composition are not yet understood. The 
usual bottom-up scenarios in which particles should be 
accelerated by astrophysical objects do not seem to pro- 
vide a convincing solution to the puzzle. The arrival 
direction of the primary particles should point to their 
sources because at such energies the intergalactic mag- 
netic field does not deviate their direction of propagation. 
However, the clustering of UHECR events observed in the 
available data is not statistically significant and therefore 
there is no evidence that they arise from point sources [3] . 
In addition, it would be necessary to overestimate several 
parameters of such sources and their acceleration regions 
in order to reach, marginally, the required energies [4] . 

The problem concerning cosmic ray sources is related 
to the necessity that they must be located in our neigh- 
borhood, since particles propagating at high energies suf- 
fer a rapid degradation of their energy. For protons or nu- 
clei as primaries, interactions with the cosmic microwave 
background should cause a loss of their energy due to 
photopion production. Such an effect should result in 
a discontinuity in the cosmic ray spectrum for energies 
above ~ 4 x 10 19 eV, the so-called Greisen-Zatsepin- 
Kuzmin (GZK) cutoff [5] . If the Auger Observatory data 
confirms that this feature is not observed, it can be shown 
[6] that protons must have travelled less than ~ 100 Mpc 
(attenuation length) in order to arrive at Earth with ener- 
gies larger than 10 19 eV. The attenuation length for pho- 
tons depends on their initial energy and it is less than 100 
Mpc for energies between 10 12 eV to 10 22 eV [7]. Since 
neutrinos have a very small cross section with nucleons 
within the Standard Model, it seems difficult that they 
could produce air showers in our atmosphere unless they 
had a yet unknown interaction [8] . 



In order to overcome such difficulties, another class 
of models have been proposed [9]. The primary parti- 
cle would not acquire kinetic energy continuously inside 
an accelerating region ("bottom-up" mechanism) as ini- 
tially thought. Instead, the highly energetic cosmic rays 
would be originated by the decay products of superheavy 
particles of cosmological origin ( "top down" mechanism) . 
For simplicity, we will consider that such particles have 
masses close to the GUT scale and would decay into 
known particles, as quarks and leptons that evolve fol- 
lowing the QCD model [10]. The quarks hadronize pro- 
ducing a small fraction of nucleons and pions that in 
their turn decay into photons, neutrinos and electrons 
and their corresponding antiparticles. Therefore, from 
the decay of such a supermassive particle it is possible 
to produce energetic photons, neutrinos and leptons, to- 
gether with a small percentage of nucleons. Depending 
on which kind of particle is the primary, different atten- 
uation lenghts can be obtained, so that one can establish 
at what minimum distances the supermassive particles 
sources should be located. 

There are different exotic candidates to play the main 
role in top-down scenarios, such as decaying topological 
defects [11] or evaporating primordial black holes [12]. 
The simplest top down models (at least from the particle 
physics point of view) involve supermassive metastable 
particles sometimes called WIMPzilla's [13]. Due to their 
colossal masses, such particles were presumably produced 
during the post-inflationary epoch and could contribute 
to a part or to the whole of the dark matter that accounts 
for about 30% of the energy density of the Universe. In 
order to explain the theoretically estimated UHECR's 
fluxes [14], such particles must be decaying now and have 
to be located in our neighborhood, which is expected, as- 
suming they are concentrated in our galactic halo. 

In this work we study the possibility of producing 
WIMPzillas in the post-inflationary process called in- 
stant preheating, suggested by Felder, Kofman and Linde 
(FKL), originally proposed as an alternative preheating 
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mechanism [15]. This process seems to be essential for 
particle production in models of quintessential inflation 
[16,17]. 

In such a scenario, scalar particles \ are non- 
perturbatively produced from the coherent oscillations 
of the inflaton 0, have their masses "boosted" due to 
their coupling to the field </>, and subsequently decay into 
supermassive metastable fermions ip. The idea of exam- 
ining stable supermassive particles in this context was 
addressed by Felder et al, but its consequences either as 
dark matter or as cosmic rays primaries were not calcu- 
lated in detail. 

More specifically, there is a relation between the den- 
sity parameter of these particles and their lifetime. If 
^-particles compose the whole of the dark matter (f2^, = 
"^"^ ~ 0.3) [181, a maximum lifetime limit will be 
found. On the other hand, a lower limit on the abun- 
dance of such particles can be obtained if lifetimes are 
constrained to be larger than the age of the Universe 
(t^ > 10 10 yr) [19]. As will be shown later, such limits 
impose severe constraints on the parameters of the FKL 
mechanism. 

This paper is organized as follows. In the next section 
we review some features of the nonperturbative processes 
more directly related to the production of massive scalar 
particles. In Section III we perform a detailed calculation 
of superheavy tp particle production, extending the pre- 
vious results in [15]. In Section IV we discuss our main 
results for produced particles considered as dark matter 
in our current Universe and present the parameter space 
for this model, which is in accordance with cosmological 
data. In the last section we present a summary of our 
main results and discuss their consequences. 



II. PRODUCTION OF X -PARTICLES 



After inflation, matter had to be (re)created, since the 
Universe became empty. In the case of chaotic inflation, 
the scenario considered here [20] , particle production may 
occur during the quasiperiodic evolution phase of the in- 
flaton field. To study such a period we assume the la- 
grangian: 



Friedmann equation: 



2 _,2 



1 



2/2 2 

9 <t> X ■ 



The inflaton field <\> produces quantum scalars bosons 
X with bare masses m x through the interaction term 
—g 2 <f> 2 X 2 For simplicity, we will limit our analiscs 
to the quadratic potential V((p) = m'i<j> 2 /2 (the simplest 
chaotic inflation model) with ss 10~ 6 Mp;, as required 
by microwave background anisotropy measurements. 

During inflation we can neglect the contribution of the 
X field and the equation of motion for the </> field is given 
by: 



■ + 3H<j) + mi 



0, 



(2) 



where H = a/a is the Hubble parameter and obeys the 



H 2 = 



8tt 
3M% 1 




m 2 4> 2 



(3) 



As far as the slow roll conditions are valid {<j> <C 3H<j>, 
4> 2 /2 < V (</>)), (j) w Mpi/3), the Universe is in an infla- 
tionary phase. Right after inflation, the <f) held oscillates 
about the minimum of its potential, with decreasing am- 
plitude due to the damping term 3H (j>, and the solution 
of the above equation becomes: 



Mpi sin (m^t) 
3 m^t 



(4) 



The <fi field may produce x-particles during this oscil- 
lating phase due to nonperturbative effects [21,22], pro- 
vided the coupling constants have appropriate values. As 
X particles are bosons, such a process may lead to an ex- 
plosive particle production through parametric resonance 
of the x ne ld [23] . To illustrate this behavior, we consider 
the quantum nature of x 

* = (2^372 / d " k ^ Xk(t)e-^ s + a k xUt)e^% 

(5) 

where a^ e a k are the creation and annihilation operators, 
respectively. The equations of motion for the fc— modes 
of the x eigenfunctions are given by: 



\',,(/)+:-!//\/,(7)4 ( ^+m 2 x +9 2 ^{t)^xk{t) =0. 

(6) 

Rewriting the above equation in terms of a more conve- 
nient variable Xu = a 3 / 2 Xk, we obtain: 



k 2 
a 2 



X k +[—+m 2 x + g 2 4> 2 )X k = 0, (7) 



where we used the fact that, for the quadratic chaotic 
potential, the inflaton coherent oscillations redshift as 
nonrelativistic matter. Note that this is an oscillator 
equation with a variable frequency 



' fc 2 
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+ m 2 +g 2 (j ) 2 {t). 



(1) The effective mass of x IS defined as 
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(8) 



(9) 



Depending on the values of the parameters, the time vari- 
ation of Uk will not be adiabatic. The loss of adiabaticity 
takes place when (f> field passes through the minimum 
of its potential, the region where u>k{t) changes more 
rapidly. In such a case, there will be an inequivalence 
between the X vacua defined before and after the pas- 
sage of the inflaton through the minimum, which can be 
interpreted as production of x particles [23]. This par- 
ticle production process has been considered mainly for 
preheating proposals, since it happens before the usual 
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perturbative reheating. Alternatively, such a coupling 
between <j> and x can be used in models with produc- 
tion of heavy metastable particles in the early Universe. 
In such a case, the produced particles may have masses 
larger than the inflaton mass [24]. In what follows, we 
work on the latter approach and consider g and m x as 
free parameters to be estimated from the available cos- 
mological data. 



III. THE FELDER-KOFMAN-LINDE 
MECHANISM 



Felder et al. [15] noticed that the preheating could be 
efficient even in a single <p field oscillation if the x par- 
ticles produced at the minimum of V(<p) were allowed 
to decay perturbatively into tp fermions when <p reached 
the maximum of its potential energy. The decay of such 
fermions into other particles and their subsequent thcr- 
malization could complete the reheating process without 
the necessity of parametric resonance. The authors also 
suggested the possibility of using such a mechanism to 
produce heavy particles in order to explain cosmic rays 
with energies above the GZK cutoff in a top-down ap- 
proach. Here, we perform explicit calculations following 
such a suggestion by assuming that tp particles are pro- 
duced non-relativisticaly. In this way the energy transfer 
process from the inflaton field is more efficient, since only 
a very "fat" x particle can decay into tp particles. The 
whole process can be schematically represented as 



9 9 i 



where g' is the coupling constant of a Yukawa interaction 
g'tptpx added to the lagrangian Eq.(l) to account for the 
interaction between tp and x- 

From the Yukawa interaction term, the decay rate of 
X bosons into tp fermions is 



T(x -> ipip) = 



g' 2 mf 
16tt 



2m,j 



cir 



(10) 



Note that the above rate is not constant, since to° de- 
fined by (9) is time-dependent, and x particles tend to 
decay at large values of (p(t). This is the most inter- 
esting feature of the model: a great amount of energy 
can be transfered from (p to the stable fermions because 
the bosons x decay when their variable masses are at a 
maximum value. 

It is possible to obtain a relation between the couplings 
g and g' if we take the interval At = tf — ti around the 
maximum value of <p evolution for which m° ff is large 
enough to allow the decay of x particles into nonrela- 
tivistic tp fermions. Denoting by and rif the number 
densities of x before and after the passage of (p through 
the local minimum of V(<p), repectively, we have that 



In 



It is convenient to work with a more intuitive time vari- 
able, i.e. the number of oscillations 



N=^ (12) 
of (f>. In such a case, the solution (4) for (p becomes 



<1>(N) 



Mpi sin (2irN) 
~~3 2irN ' 



(13) 



This is a good parametrization from N = 0.25 and later 
if we fix N = 0.5 to be the instant when the inflaton field 
crosses the minimum of its potential for the first time. We 
will consider that the maximum momentum of ip particle 
to be its mass, that is, p^ < m^,. The maximum avail- 
able energy for the creation of a ip pair is <?<E>, where $ 
is the amplitude of (p(N). This determines the largest 
interval AN around the maximum amplitude of the <p 
field, N w 0.72 ($ w 0.07M P; ), for which mf w g<P is 
large enough to allow a pair creation. Such an interval 
is found to be A./V = 0.25. Additionally, we require that 
approximately 90% of the x particles decay during this 
interval. From Eq. (11) we obtain, 

ln(^)<',^ r^<p(N)dN, (14) 

with <p(N) given by (13). Solving numerically the above 
integral, we obtain that 



g' 2 9 > 3 x 10" 



(15) 



T(t)dt. 



(11) 



This upper limit is almost one order of magnitude larger 
than the Felder et al estimate [15] (g' 2 g w 5 x 10~ 4 ), 
since they did not consider the phase space factor in the 
expression for the decay rate (10). It is important to 
notice that the only arbitrariness in our assumptions is 
the fraction of remaining x particles (90%), but for rea- 
sonable choices (say nf/rii between 1/2 and 1/100) the 
constraints on g and g' do not vary appreciably. All the 
other constraints are consequences of the assumption of 
a maximum momentum p^, < m^, so that tp may be con- 
sidered as non-relativistic. Naturally, the limits on g' 2 g 
would be even tighter if we had considered values for 
larger momenta. Therefore, we are being conservative in 
our estimates. We also verified that although one should 
assume the time dependence of the decay rate on the 
above calculations, it does not bring any important dif- 
ference if compared to the estimate found in [15] where a 
constant decay rate, T, was used. This seems reasonable, 
since the integration interval for N is taken to be around 
the maximum of (p(N) where the sine function does not 
vary significantly. Such care would be necessary if we 
were studying the production of relativistic particles. 

Independently of the calculation details, our main goal 
in this work is to find the largest possible region of the 
parameter space of the FKL mechanism that is phe- 
nomenologically allowed by the available data. Assuming 
that the tp fermions are the metastable massive particles 
that we are looking for, we need to evaluate their present 
abundance supposing that they contribute to the energy 
density of the dark matter today. Since for each decay- 
ing x a tptp pair is created, we have to find the number 
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density n x of the \ bosons for each oscillation of the in- 
flaton. For the first oscillation, n x can be calculated from 
the solution of Eq. (7) about the minimum of V{4>), so 
that the total number of x-particles [15] is 



1 



I 



d k nk = 



8tt 3 



■ exp - 



g\<Po\ 



(16) 



and <j>o = rri0$ is the field velocity near the minimum 
of the potential. The model describes the above produc- 
tion of x particles, the boosting of their masses and their 
subsequent decay into tp WIMPzillas with masses 



w 5$ w 0.07gMpi. 



(17) 



In order to verify that most of the tp particles will be 
produced in the first oscillation of (p, we compare the 
number density of produced x particles in the second 

(2) 

passage by the minimum of the <p potential, n x , to the 



first one, n x ^ . By taking into account the dilution of n x 
due to the Universe expansion between the consecutive 
passages, it is found that (for details see the Appendix) 



(i) 



n (2) 



= V2 



exp 



-6tt 2 x 10 6 



(18) 



Therefore, the exponential term could provide the desired 
suppression between the two first passages. In fact, it can 
be shown that the ratio between two consecutive passages 
tends to a constant prefactor multiplying the exponen- 
tial, and, as long as the suppression is assured, particle 
production will be negligible for all subsequent oscilla- 
tions. We will see in the next section that such suppres- 
sion is verified, since for typical values of the parameters, 



say, 



g = 1CT 2 and m x = 10~ 4 M P; , n 



(2) 



10" 



26 n (D 



and 



hence it is reasonable to consider only the first passage 
in our calculations. 



IV. tp ABUNDANCE AND ULTRA HIGH 
ENERGY COSMIC RAYS 

If the tp particles are the superheavy relics that decay 
into the observed UHECR's, we can use the presumed 
cosmic ray flux associated to their decays [2] in order to 
estimate limits on the parameters of the model discussed 
above. 

By accounting for the dilution of tp particles since their 
production until today, one may find the associate den- 
sity parameter O^(fo) = pi/j(to)/pc- In order to obtain an 
expression for this parameter, we have to consider three 
different moments in the history of the Universe: the 
production of ip particles (t p ), the end of the reheating 
period (t r h) and today (t ). We can then write: 



Pn{to) 



PR{trh) 



(19) 



where we assumed thermal equilibrum for the relativistic 
components and the fact that p^p(to)ao 3 = Pip(t r h)arh 3 ■ 
On the other hand we can suppose that the Universe will 



be reheated from an instant convertion of the remanes- 
cent inflaton density energy into relativistic particles, so 
that 



Pj){trh) P^jUh) _ Pip (t p ) 
PFl(t r h) P4>{trh) P<l>{tp) 



(20) 



where the last equality is obtained by considering that 
the tp and <p field coherent oscillations redshift as non 
relativistic matter. Note that we are supposing that only 
the tp particle production in the first passage is important 
(see expression (18)). 

Finally we can substitute p<p(tp) = (iM Pl /%TT)H 2 into 
(19) so that 



n^{T )h 2 



87r m^riif, 



(21) 



where Tq w 2.7-K" is the present CMB tempera- 
ture, h = ^/(lOOkm.s-i.Mpc- 1 ) and n R (T )h 2 w 
4.3 x 10~ 5 is the current radiation density parameter. 
The other parameters are H 2 = (8n/3M Pl )(p < p p ) — 



{%n/ZM Pl ){m^/2) « (47T/3) x 10 



-14 M 2 

ivi pl 



and, given 



that n x — and for nonrelativistic fermions m c x 



off 

m x n 



2m,/,, we find that m^,n^, 
by equation (16). Choosing T r h 
equation can be rewritten as 



where 



n x is given 



X' 

10 9 GeV, the above 



n^{T )h 2 w 4.72 x 10 18 



(. 



T, 



rh 



\ 10 9 GcV 



'2.7K 



.3/2 



10 



~ 2 g 2 exp 



7T X 10 7 



(22) 



If tp WIMPzillas arc required to explain UHECR's, it 
is possible to obtain a relation between their abundance 
(fi^,), their masses (m^), the lifetime of tp (r^) and the 
UHECR's fluxes (for details see references [2,19]). As 
we mentioned earlier, in order to obtain such relations 
it is necessary to adopt a model in which particles are 
produced from the tp decay. This topic is an issue by 
itself, and we will follow the usual assumption of extrap- 
olating QCD mechanisms valid in the the quark-hadron 
fragmentation process for the higher energies considered 
here. As we will show later on, the use of a specific value 
for the flux will not alter the region of the allowed pa- 
rameter space appreciably. Such a choice implies that 
photons dominate the primary spectrum by a factor of 
w 6 over protons [2]. This means that, if we consider the 
observable ultra high energy cosmic ray fluxes as due to 
extremely energetic photons resulting from the decay of 
tp particles, we can write: 



t> = 3.16 x 10 i8 / {Q^h 2 ) 



M 



pi 



1/2 



(23) 



where / measures the clustering of tp particles inside the 
galactic halo. It is taken as 1 for a uniform distribution 
of superheavy WIMP's, but can be considered as approx- 
imately 10 3 if they are concentrated in galactic halos (see 
[2]). We will assume the latter case in what follows, i.e. 
/ ~ 10 3 , so that following limits on the parameters of 
the model can be established: 
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• if they constitute the whole of the dark matter (for 
O/,~0.3, ft ~ 0.7) [18]) 

^{T )h 2 



0.14: 



(24) 



• if they have a lifetime of the order of the age of the 
Universe 



7V, w 10 lu years -» O/,(T )ft 2 w 3.16 x 10 



-12 



1/2 



(25) 

Therefore, the V abundance obeys the limits 3.16 x 
10- 12 (.g$/M P; ) 1/2 < O^ft 2 < 0.14. This imposes con- 
straints on the x bare mass according to (22): 



3.2 x 10 6 g In (^4.72 x 10 19 g 3 / 2 (j£jh 
3.2 x 10 6 g In (u49 x 10 30 g 3 ' 2 



10 



V 



f) 



10 



V 



1/4 



(26) 

Such conditions define a parameter space which is conve- 
nient to assess the phenomenological viability of the FKL 
mechanism as a generator of ip non relativistic WIMPzil- 
las. This is possible because the number density of \ 
produced in the post-inflationary era depends on g and 
m x and is related to Q^,, so that limits on today re- 
strict the possible values of such parameters. The analy- 
ses is summarized in Figure 1, which can be understood 
as follows. The upper and lower limits of the gray area 
come from the substitution of Eqs. (24) and (25) in (22). 
We also have imposed unitarity constraints on Eq. (15) 
(g and g' < 1), so that 3 x 10~ 3 < g < 1 (which com- 
prises the upper right triangle in the figure) and obtained 
the limits to the left (the minimum value for g given the 
maximum possible value for g') and to the right (the 
maximum value for g) in the allowed (dark gray) area in 
the figure. For the sake of comparison, we also included 
the same analysis for the original FKL result, which cor- 
responds to the right triangle limited by the dot-dashed 
line and to the widened area including the light gray re- 
gion. 




Ml 
O 



FIG. 1. The allowed parameter space (m x ,g,g') for the 
production of WIMPzillas in the FKL mechanism. For m x 
given in Planck mass units (m Xrl ) , the dark gray area stands 
for the allowed values of the parameters between h — 0.14 
(dashed line) and fl^h 2 = 3.16 x KT 12 (g<S>/M Pl ) 1/2 (solid 
line). Unitarity constraints on g' and g provide limits to the 
left and to the right in the allowed area in the figure. The 
application of these limits on the original FKL results [15], 
g' 2 g ~ 5 x 1(U 4 , increases the allowed area by adding the 
light gray region. 



We see that the allowed region in the parameter space 
is rather constrained. Particularly, given the valid range 
for g, we find very high masses for ip. From Eq. (17) 
the minimum ip mass that can be obtained through 
this model is ~ 10 15 GeV (the upper limit being 
~ 10 18 GeV. This happens because the exponential 
suppression of the number density of \ particles created 
after the first <p field half oscillation, necessary to avoid 
parametric resonance (see Eq. (18)), is also present in 
the expression for Cl^ (see Eq. (22)). Although the upper 
cosmological limit is quite strong given the most recent 
measurements, one may consider i\ir> <X weak constraint 
the second astrophysical cosmic ray limit, as it is model 
dependent. However, even if we consider other classes 
of models to establish new lower limits on fi^ the above 
results would not change significantly, because the pa- 
rameters of the model are not very sensitive to variations 
on O^,. Let us assume, for the sake of a comparison, 
an hypothetical value for the present ip abundance, say 
Q^h 2 <~ 10~ 20 , so that the allowed region of the Fig. 1 
is enlarged. The resulting parameter space is depicted in 
Figure 2. 




o 



FIG. 2. The allowed parameter space of the previous fig- 
ure is not substantially altered even if it is enlarged by an 
exagerated and hypothetical lower limit for the abundance of 
ip (fl^h 2 — 1(U 20 , solid line). All the other lines have the 
same meanings as in Fig. 1 

For the same reason, wide variations of the reheating 
temperature in Eq. (21) will not change the picture. For 
example, one would have to consider reheating tempera- 
tures about 20 orders of magnitude higher than the one 
assumed here to obtain a shift of only one order of mag- 
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nitude on the allowed range of m Xpl . Therefore, the re- 
sulting allowed area of the parameter space is relatively 
independent of particle physics details of the reheating 
and of the hypothetical top-down decay of tp particles in 
UHECR's. 



V. CONCLUSION 



We studied how to generate supcrmassive fermions 
that can explain the UHECR's in the context of a par- 
ticle production mechanism suggested in Ref. [15]. We 
obtained the parameter space for which such a mecha- 
nism can take place and concluded that some fine tuning 
of parameters seems to be necessary. Additionally, the 
lower limit on the tp mass obtained, m$ > 10 15 GeV, is 
rather robust. A typical signature of this model would be 
an unforeseen rise in the flux at the highest energies end 
of the cosmic ray spectrum, which could be observed by 
the next generation of experiments, like the Pierre Auger 
observatory. The cosmic ray spectrum must have a cut- 
off which is associated to the maximum energy possible 
to UHECR's and is independent of the GZK feature. If 
such a cutoff happen to be below 10 15 GeV, this simplest 
version of the FKL mechanism should be discarded. On 
the other hand, it is this mechanism that can provide 
masses of such magnitude more naturally than any other 
top-down versions, so that if it is at all possible to mea- 
sure such high energy cosmic rays and no cutoff in the 
UHECR's spectrum is observed by the next generation 
experiments, this model can become an attractive can- 
didate. It is also interesting that, despite having per- 
haps too many free parameters, this model is rather con- 
strained, and such a result is relatively insensitive to wide 
variations of the relevant cosmological parameters. 

It is important to emphasize that we studied the pro- 
duction of non relativistic tp particles only, and the sce- 
nario can be made more complex by considering the pro- 
duction of tp particles that are relativistic at the pre- 
heating time but becomes non-relativistic along the Uni- 
verse evolution. In this case, the energy transfer from 
the inflaton field to other fields may be not very effi- 
cient and it would be necessary to consider the dilu- 
tion/concentration of tp particles throughout the several 
phases that happened since inflation (coherent oscilla- 
tions phase, radiation domination and matter domina- 
tion) and the allowed parameter space may be widened. 
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APPENDIX A: THE EXPONENTIAL 
SUPPRESSION 

From (13), we can write the time derivative of <p(N) 
with respect to N: 



<f>'(N) 



Mpi 2ir cos (2ttN) M Pl sin (2nN) 
~ 2k~N 



(Al) 



3 2wN 3 2irN 2 

By labelling each time that <p passes through the mini- 
mum of its potential as Nj we have that 



Mpi cos (2wNj) 
~~3 Wj ' 



(A2) 



that is, 



1 passage : 
Nj. = 1/2 

2 nd passage 
jV 2 = 2/2 

3 rd passage : 
N 3 = 3/2 

j th passage : 
Nj = 3/2 



= m 



PI 



\M Pl 



Ws\ = 



1^1 = 



lM Pl 



Mpi cos (Jn) 
3 j/2 



From the general expression above, it is possible to ob- 
tain the ratio between the particle number density of x 
particles produced in two consecutive bursts of \ gener- 
ation so that we can evaluate the amount of suppression 
for each passage j. Writing the expression (16) for the 
total number density of \ particles produced in each pas- 
sage in terms of the new definitions: 



8tt 3 



3/2 /|<//. k 



2tt 



3/2 



cxp 



Mpi m,p \ 
3 2irj) 



3/2 



cxp 



g\(p' j m 4> /2Tr\ ^ 
67r 2 j m 2 / 



M. 



(A3) 



Since the a(t) ex i 2 / 3 along the coherent oscilations phase, 
we must compare the number density of particles pro- 
duced at tj until tj+i by taking into account the dilution 
of rij at the (j + l)th passage: 



n«(t 3+1 )=n«(i^ 1 



(A4) 



However, from (12), tj/tj + i 
considered stands: 



j/j + 1, so the ratio to be 



(tj+i) a + 1 



,0) 



(*i+i) 



(A5) 



By using the expression (A3) in the above relation, we 
finally find: 



Si- 



.7+1 J 



n ( i\t J+1 ) 



3+1 

3 



1/2 



6^ 2 



mr 



exp 



(A6) 

Since \ particles decay into two tp's, the ratio above is 
also valid for tp particles. Thus, we conclude that the 



6 



exponential suppression, that has to be addressed in or- 
der to obey cosmological limits on the WIMPzillas abun- 
dance (see section IV), guarantees that we can consider 
in a good approximation only the x particles produced 
in the first passage. 
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